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Strong and fast "bang-bang" (BB) pulses have been recently proposed as a means for reducing 
decoherence in a quantum system. So far theoretical analysis of the BB technique relied on model 
Hamiltonians. Here we introduce a method for empirically determining the set of required BB 
pulses, that relies on quantum process tomography. In this manner an experimenter may tailor his 
or her BB pulses to the quantum system at hand, without having to assume a model Hamiltonian. 
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I. INTRODUCTION 



Quantum computers hold great promise in solving cer- 
tain computational problems faster than their classical 
counterparts, but they are notoriously susceptible to de- 
coherence (deviations from unitary dynamics) and uni- 
tary errors, the combination of which we refer to as 
"noise" . The effect of decoherence is to induce computa- 
tional errors that destroy the quantum speed-up: a deco- 
hered quantum computer can be efficiently simulated by 
a classical computer Q. Hence the ultimate success of 
quantum information processing depends on the ability 
to implement error correction or avoidance techniques. 
To this end, a variety of quantum error correcting codes 
(QECC) and other methods have been designed. These 
methods all share an important feature: they are de- 
signed to deal with specific models of errors, as embod- 
ied in an assumed system-bath interaction Hamiltonian. 
The class of active (e.g., stabilizer) QECC § | |, |, 
for example, is designed to correct independent errors 
resulting from (up to) some fixed number, t, of system- 
bath many-body interactions; the class of passive QECC 
(decoherence-free subspaces) works optimally under the 
assumption of collective (i.e., fully correlated) decoher- 
ence H [?], ||, [| or assumes multiple-qubit errors ; 
dynamical symmetrization methods assume baths with 
relatively long correlation times and weak system-bath 
coupling, so that decoherence may be suppressed using 
fast and strong "bang-bang" (BB) pulses, introduced in 
t!|,_andjurther _developed in [p| ,Jlj pj [l6[ ^ |J, 




H|, |2g £3|, H H § In spite 

of this impressive arsenal of methods there is a funda- 
mental problem in the model-specific approach in terms 
of its applications to experimental quantum information 
processing. The problem is that in real world applica- 
tions, decoherence is often a combined effect, which arises 
from a variety of sources, and does not correspond to 
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one particular model. It is often very difficult to iden- 
tify and isolate the various sources. The result is that 
the model-dependent approach for overcoming decoher- 
ence breaks down when applied to realistic systems, since 
it inevitably fails to capture all sources. In addition, 
current methods tend to ignore the experimental con- 
straints imposed by the finiteness of resources, such as 
the scarcity of qubits in present-day implementations of 
quantum computers (presently, fewer than 10 qubits). 
Of course, this criticism in no way diminishes the impor- 
tance of the model-specific approach: it is through that 
approach that ground-breaking new results have been ob- 
tained which establish the in-principle possibility of over- 
coming decoherence. In particular, this work has led to 
the observation that fault tolerant quantum computation 
is possible in the independent errors model provided the 
fidelity of gate operations is above a certain threshold 

@ flung. 

We focus here on the BB method and consider a 
paradigm that is the reverse of the model-dependent ap- 
proach to decoherence-reduction: Instead of assuming a 
specific model of decoherence, designing a correspond- 
ing QECC, and then looking for a system that might 
be described to a good approximation by that model, 
we propose to tailor a set of BB pulses to a system, 
from experimentally measured decoherence data. We call 
this approach, which we introduced first in |26y , "Em- 
pirical Bang-Bang" . Empirical BB is a phcnomcnolog- 
ical approach which forsakes a microscopic understand- 
ing of the underlying decoherence processes in favor of 
a direct attack on the combined effect of all sources 
of decoherence at once. The procedure can be itera- 
tively optimized using a closed-loop learning algorithm 
PH pfj| , p7| pi |39| , f4(i| , |4l[ | . In this manner one may take 
into account practical constraints imposed by the specific 
physical and experimental realization. 

That empirical BB is feasible in principle follows from 
two key facts: (i) It is possible to experimentally mea- 
sure the superoperator (i.e., the map that propagates the 
density matrix) characterizing the noise in a particular 
system by using Quantum Process Tomography (QPT); 
(ii) As we show here, given knowledge of the superoper- 
ator it is possible to design a BB procedure. Thus an 
experiment can, in principle, provide all the information 
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needed to design an optimized set of BB pulses. 

This article is arranged as follows. In section II we 
review the basic background to quantum process tomog- 
raphy, and the theory for decoupling by symmetrization. 
We then present, in section a derivation and discus- 
sion of several formulas for determining the set of decou- 
pling operations. These results are illustrated in section 
IV with a few examples. We then indicate in section |v| 
how the empirically determined set of BB pulses can be 
optimized using a learning loop algorithm. 



II. REVIEW 

In this section we review the important components 
of the empirical determination of bang-bang operations. 
These include quantum process tomography (QPT), the 
theory of dynamical decoupling operations for a given 
Hamiltonian, and its geometrical interpretation. Readers 
familiar with these concepts can choose to skip ahead 
to section III, although the notation introduced in this 



section will be used in the remainder of the paper. 



A. Quantum Process Tomography 

The dynamics of an open quantum system coupled to 
a bath is formally obtained from the time-ordered evolu- 
tion 



U(t) = TexpH / H(t')dt') 



under the combined system-bath Hamiltonian 
H = H S ®I B +I S ®H B +H SB 



(1) 



(2) 



where / is the identity operator, H$ is the Hamiltonian 
for the system alone, H B is the Hamiltonian for the bath 
alone, Hsb is the system-bath interaction Hamiltonian, 
and the 5* 7 and B 1 are operators on the system and the 
bath respectively. Tracing over the bath degrees of free- 
dom in order to obtain the time-evolved system density 
matrix: 



P (t) = Tc B [U(t) (p(0) ® MO)) U*(t)], 



(3) 



where p(0) is the initial density matrix of the (open) sys- 
tem, p B (0) is the initial density matrix of the bath. It can 
be shown that this agrees with the most general quantum 
evolution consistent with the condition of complete pos- 
itivity, known as the Kraus operator sum representation 
(OSR) HHg: 

£ t ( P (0)) = P(t) 

= E X«At)K a p(0)Kl- (4) 



The Kraus operators can be related to Eq. (H) through 



(5) 



^V(t) = y/\ v {n\U(t)\v) 



where \v) 1 \p) are eigenvectors of the the initial bath 
density matrix: /5b (0) = J2v @- Since 

Tr[p(i)] = 1, they satisfy the normalization condition: 
YtuMivApv = I s - The matrix 



fit/ 



is a time-dependent, Hermitian coefficient matrix defined 
by a transformation of the Kraus operators to a fixed (i.e., 
time-independent) operator basis K a : 



A prescription for determining the superoperator £ t from 
experimental data (QPT) was given in a number of recent 
papers |^], fl8| , and has very recently been applied in 
NMR experiments p9[ . In this paper we will take QPT to 
mean the determination of the coefficient matrix Xa@{t)> 
with respect to a given (experimentally convenient) choice 
of fixed basis operators K a . Formally, the problem is 
to invert the ^-matrix from experimental data. Since 
X is time-dependent it is clear that one can in practice 
only sample it. If the decoherence process is Markovian 
then it suffices to obtain the time-independent coefficient 
matrix A that appears in the Lindblad equation |^0|, ^l) . 
However, even this is a formidable problem: if the density 
matrix has dimensions N x N (where for n qubits N — 
2") then a simple counting argument shows that there are 
at most iV 4 — V 2 independent real parameters in A and 
the same number, but time-dependent, in x- Even for 
one qubit this amounts to 12 different parameters that 
may have to be measured to completely characterize the 
decoherence process. Fortunately, it is well known that 
in practice as few as 2 parameters may suffice, as is the 
case with the T\ and T2 relaxation times in NMR fl52[ . 

The general idea behind QPT is to characterize the su- 
peroperator action on a complete basis set. To see this, 
let the N 2 matrices pj be a basis for the density ma- 
trix p. For example, pj could be the set of pure states 
\j)(j'\, which are then fed into the decoherence process as 
inputs: £(pj) = J2k ^jkPj- Using quantum state tomog- 
raphy f33f , one can experimentally determine Xjk, which 
fully specifies the superoperator £ , since it is now possible 
to find the x-matrix: Define £ by K a pjK^ = 'Y2 k £,^ Pk, 
where K a are the fixed basis Kraus operators. Then 
one can show that J2 a p ^k Xap = Xjk Hi- This can 
be thought of as a matrix equation for the vector x and 
it can be solved by computing the inverse of the ma- 
trix £. Thus, by measuring A and by giving £ through 
a choice of the fixed operator basis K a , finding the x~ 
matrix has been transformed into a linear algebra prob- 
lem. In practice, we note that it may often be difficult to 
prepare the full basis set pj. An interesting alternative, 
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using entangled input states, was recently proposed in 
Q. A method that circumvents tomography altogether 
(but is less general), using quantum network ideas, was 
described in |55l. 



B. Decoupling by Symmetrization 

The process of decoupling by symmetrization counter- 
acts noise by applying sequences of frequent and strong 
pulses. The time scales are crucial: one needs to per- 
form a complete cycle of symmetrization operations in 
a time shorter than the inverse of the_ high-frequency 
cutoff of the bath spectral density 
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An 



elegant group-theoretical treatment shows that the ap- 
plied pulses are unitary transformations forming a finite- 
dimensional group, and the application of a series of 
pulses amounts to an average (symmetrization) over this 
group jfl| [H| [IT], pi p^|. A geometrical interpreta- 
tion, reviewed below, can offer further insight p4j. The 
method can also be used to perform "environment engi- 
neering" , in order to prepare the conditions that allow 
for DFSs (l^, [ll| |?J, as well as in order to eliminate 
leakage errors that couple encoded states with states out 
of a DFS ITF ' 



28 



We briefly review this theory. 
A set of symmetrization operations is chosen such that 
they form a discrete subgroup of the full unitary group 
of operations on the Hilbert space of the system. Denote 
this group G and its elements gj, j = 0,1,..., \G\ — 1, 
where \Q\ is the order of the group. The cycle time is 
T c = \Q\At, where \Q\ is the number of symmetrization 
operations, and At is the time that the system evolves 
freely between operations under Uq. The symmetrized 
evolution is given by 



ISI-l 



U(T C ) = J] g}Uo(At) gj 



where the evolution under H$ b + Hb has been neglected 
during pulse application, i.e., during the action of the 
group elements gj. H e s denotes the resulting effec- 
tive Hamiltonian. Since the approximation requires very 
strong, short pulses to be implemented in a sequence, 
they have been termed bang-bang (BB) operations (we 
will use decoupling, symmetrization, and BB operations 
interchangeably). In this (BB) limit 

ISM 

h ~ H eS = — g) H 9 3 = n e (H), (6) 

11 3=0 

where H c g is the desired Hamiltonian (without noise). 
The map is the projector into the centralizer, Z(Q), 
defined as 

z(g) = {x\ [i !Si ] = o,v ffj eg}. 

It is clear that Hg commutes with all gj so that, if 
our group is generated by {/, Hs, £7}, the evolution will 



proceed without the operators S y affecting the system 
since the error operators will commute with the effec- 
tive Hamiltonian. The control algebra is the algebra gen- 
erated by the set {gj}. Even if the symmetrization is 
performed under less than ideal conditions, BB can still 
reduce the noise in the system 

The main advantage offered by dynamical decoupling 
is that it does not require extra qubits. This is a very 
attractive feature compared to both active and passive 
error-correction, one that may make dynamical decou- 
pling a method of choice for small-scale quantum com- 
puter implementations, provided its stringent time-scale 
requirements can be met. 



C. Geometry of the Decoupling Method 

In preparation for the remainder of the paper, and as 
an intuitive aid, we briefly review the geometric descrip- 
tion of BB controls developed in p^] . Let us explicitly 
introduce N = n 2 — 1 traceless, Hcrmitian generators 
{^i}iLi °f SU{n). These generators are closed under 
commutation and span the space of traceless Hcrmitian 
matrices. For SU(2), the Pauli matrices are commonly 
used; for SU(3), the Gell-Mann matrices, and for higher 
dimensions, one may use a direct generalization of the 
Gell-Mann matrices. For dimensions that are a power of 
two (and quantum computing) it is often convenient to 
use the Pauli group (tensor products of Pauli matrices). 
The {Ai} satisfy trace-orthogonality, 



Tr(AiAj) = MS 



13 > 



(7) 



where M is a normalization constant (often taken to be 
2 for Lie algebras or n for n x n matrices). Expanding 
the system operators in terms of the {A^} yields: 



where the expansion coefficients are 

Using this, H$b can be written as as follows: 

N 

h S b = ff 7 ® B~ t = a ll \ l ® b 1 

7 7 i— 1 



(8) 



(9) 



^(a 7 • A) Cg) .B7 



(10) 



Here a 7 and A are vectors of length TV. In this represen- 
tation, used extensively in ]5(|, an n x n Hamiltonian, 
if, is a vector with coordinates a 7 for each error 7 in an 
A^-dimensional vector space spanned by the {A^} as basis 
vectors, with ordinary vector addition and scalar multi- 
plication. The open system evolution is thus described 
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by a vector (or vector field) in the space of possible evo- 
lutions. 

Now, as is well-known, there is a homomorphic map- 
ping between the Lie groups SU{2) and 50(3) J57|. This 
mapping is generalized as follows for SU(n) and a sub- 
group of the rotation group SO(N): 



N 



(ii) 



where the matrix R^ £ SO(N), the adjoint representa- 
tion of SU{n). 

The BB operation [Eq. (0)] may now be viewed as a 
weighted sum of rotations of the (adjoint) vectors a 7 . To 
see this, first let 



R™cL 



(12) 



This represents the rotation by R^ of the coordinate 
vector a-y. Next average over all rotations: 



i \e\-i 

— V d {k) 
\Q\ ^ 7 ' 

11 fe=0 



(13) 



Finally, note that the effective Hamiltonian, after the BB 
operations, can be rewritten as: 



161-1 

|y| fc=0 7 



(14) 



Eq. (O) [compare to Eq. jlC|)1 is the desired geometric 
representation of BB operations. Their effect is to simply 
transform, for each error 7, the coordinates a 7 to a'. It is 
simplest to interpret this in the case of storage, where we 
seek BB operations such that H c g = 0. Since the errors 
can be decomposed in the linearly independent basis set 
indexed by 7, each term a 7 • A must vanish separately. 
Furthermore, since the Xi are independent this can only 
be satisfied if a' — for each 7. This means that 



(15) 



i.e., the sum of all rotations applied to the original coor- 
dinate vector a 7 must vanish. 

Similarly, to obtain a modified evolution correspond- 
ing to a target Hamiltonian i?* ff = X) 7 (^7 ' ■*) ® ^11 we 
require the weighted sum of rotations applied to the orig- 
inal coordinate vector to be equal to the corresponding 
target coordinate vector SL. I.e., for iJ c ff ^ 0, the fol- 
lowing condition should be satisfied to obtain the desired 
evolution: 

^ - < (16) 

This may require a combination of switching strategies 
for the BB pulses pi. 



It should be noted that the geometrical picture is an 
explicit representation of a subset of the group algebra 
CQ using the set of traceless Hcrmitian matrices and the 
identity as the basis. When the coefficients of the adjoint 
vector are real, the resulting matrix H e g is Hermitian. 
When they are complex, the resulting matrix is not Hcr- 
mitian and the evolution is not unitary, but may still be 
treated empirically, as we show below. 

We now turn to showing how to find the BB pulses 
directly from experimental data, i.e., given a QPT mea- 
surement of the ^-matrix. 



III. DETERMINATION OF SYMMETRIZATION 
OPERATORS 

Since the BB method operates at extremely fast time- 
scales it is useful to consider a short-time expansion of the 
OSR evolution equation dJ). To do so we follow jl5[ |58| , 
where it was shown how the OSR can be rewritten to 
resemble the Lindblad equation 50 0. Thus, the OSR 
can be rewritten as 



P(t) = P (0)--[S(t),p(0)] 

N 



a, 0=1 



\ £ xaAt) ([K a , P {Q)Kl] + [jr Q p(o), jrj]) , 

(17) 



where S(t) is the Hermitian operator defined by 



S(t) 



ih 



N 



mi] 



(18) 



a=l 



Note the similarity of Eq. ( O ) to the Lindblad equation 
p(j| . |5~i[ . Indeed, the Lindblad Markovian semigroup mas- 
ter equation can be derived from Eq. dl7| ) via a coarse- 
graining procedure |45| , |58| , which replaces the time- 
dependent Xa,p matrix elements with their time-averages 
over an interval that is longer than the bath correlation- 
time, and thus longer than the BB time-scale. An im- 
portant outcome of this procedure is that the coarse- 
grained S(t) can be interpreted as a system Hamiltonian 
Hs plus a Lamb shift correction |45|, Q . While still ex- 
act, Eq. ( Jl7[ ) is more amenable to a short-time expansion 
than the original (equivalent) form of the OSR, Eq. (^). 

Note that the "fixed-basis" {K a }a=i ^ s completely 
analogous to the Hermitian generators { Ai}^ 1 of SU(n) 
used in the geometric picture of section |II Oj . Thus, as- 
suming a Hermitian basis {K a } we can rewrite Eq. Jl8| ) 
as 

M 

S(t) = thJ2Mx a ,o(t)}K a = ihlm(x) ■ K, 

which can be interpreted as giving the "Hamiltonian" 
S(t) as a vector with coordinates Im[xa,o(£)]} m a space 
with basis vectors {K a }. 
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Next we give a general method for determining BB 
controls from empirical data, specialize the applicability 
of this method somewhat, and then treat storage, single 
qubit operations and computation. 



A. Empirical Bang-Bang Condition 

Before going into a detailed and more careful analysis, 
we first present a "rough" version of the empirical BB 
condition. We note two key facts: (i) the BB method will 
operate only to undo the undesired evolution due to S(t); 
(ii) from sections II A and 



II B 



we find that, under the 



action of BB controls, the {K a } transform as 



|y| k 0=1 



(W) 



Thus, given the considerations above concerning the ef- 
fect of BB pulses and their geometrical interpretation, we 
can express the BB-modified open system evolution as 

S = ihlm(x) ■ K 

|y| k a/3 

= ifilm(x) K = S, (20) 

where the new, BB-modified "Hamiltonian" S is de- 
scribed by the new, rotated coordinate vector 



Im(|)=Im(x)~^i? (fe) . 



(21) 



Now, let the ideal, or desired "Hamiltonian" be described 
by the coordinate vector \ w , i.e., 



S w = ihlxsx(x w ) ■ K. 



(22) 



For storage this would correspond to the null vector, but 
not for computation. The goal of the empirical BB pro- 
cedure is to find rotation matrices such that the 
difference 

S - S w = ih [lm(|) - Im(x™)] • K = 0, (23) 

or more generally, is minimal. This has the simple ge- 
ometric interpretation of minimization of the distance 
between the BB-modified vector Im(x) and the desired 
vector lm(x w ). 

The input data is Im(%) (the output of the QPT mea- 
surement), Xw (the desired Hamiltonian), K (the opera- 
tor basis, with respect to which Xw and x are defined). 



This data specifies a solution to Eqs. (pl|), ( p3| ) in terms 
of the rotation matrices R^ . T his sol ution is not unique; 
see, e.g., the example in section IV A. 

When the RS k ' are found, the BB pulses can be cal- 
culated from the transformation connecting the adjoint 
representation to its unitary group. 

Thus Eqs. (pH),(p^ can be viewed as the essence of the 
empirical BB procedure. From here on we flesh out this 
first main result. 



B. Qubit Noise 



The development in section [II A was cavalier in its 
treatment of the indices a, [3 of the fixed operator basis 
{Ka}. To be more precise, consider a quantum regis- 
ter of N qubits. We will derive a short time expansion 
of Eq. Jl7 ) under the assumption that the system-bath 
interaction is linear in the system operators: 



N 



rr(l) 
SB 



Bi 



(24) 



where <7i = (af , a\ , of) is the vector of Pauli matrices 

acting on the i th qubit, and Bi is a corresponding vector 
of bath operators. This assumption will be relaxed below 
(Section |III E| ) and, as should be clear from section III A 



is not essential for our approach, but will make the cal- 
culations below more transparent. A Taylor expansion 
of the evolution operator U(t) [Eq. (0)] then reveals that 
as time increases, higher and higher tensor powers of the 
Pauli matrices act on the qubits: 



U(t) = I - itH 



(i) 

SB 



2! 



(25) 



where for simplicity we have assumed a time- independent 
Hamiltonian and set H$ = Hb = 0. The 0(t) term 
involves only single Pauli matrices, but the 0(t 2 ) terms 
and higher involve tensor products of Pauli matrices. To 
capture this in terms of the OSR we expand the fixed 
basis operators K a as 



(26) 



where, for the i qubit, erf, a = 0, 1, 2,3 corresponds to 
Ii,af , a\ ,a* respectively. The subscript on Ks n denotes 
a vector a n — (ai, ...,a_/v) with n non-zero entries. I.e., 
Kg n acts non-trivially on n qubits. (We also use a for a 
vector of arbitrary index ) Note that we have omitted the 
subscript i on a in Eq. ( J2q ) in order to reduce the index 
clutter. There exist M — 4 N different operators 
with Ks = I ® • • • 8) / being the identity on the space of 
all qubits. Here we have chosen the if 's to be Hermitian, 
and trace orthogonal: 



2 N d. n 



(27) 



6 



Hence they are a valid basis for all 2 N x 2 N matrices. 

Corresponding to this expansion of the fixed-basis op- 
erators, we can rewrite the OSR, Eq. (0) , more explicitly 
as 



N 



Pit) 



m,n— ft 



a m ,}3n 



Xa ,A(*)p(°) 



JY 



E E X Sm< o (t)K Sm p(O) + X^^WKO)^ 



m=l <S„ 
iV 



+ E E ^A^ft.^. (28) 

Thus terms that contain only single Pauli matrices but 
not tensor products of Pauli matrices can only come from 
the second sum E TO= i Es )j with m = 1. Comparing 
to Eqs. (|l7|),(p^|) it is clear that this sum is responsible 
for (part of) the "Hamiltonian" S(t), whereas the third 
sum generates the Lindblad-like term in Eq. (|l7|). (This 
can also be verified directly by repeating the derivation 
in ^8) using the K$ n .) Hence to first order in t we 
find 

p(t)-p(0) «-i[S(f),p(0)] 

Si 

= -t[Y^lm( XSu o)K Sl , P (0)} (29) 

Si 

where in the last line we used the hermiticity of the K 
operators. The term Y] Sl is a sum over all elements 
of the Pauli group with one non-identity element in the 
tensor product. By comparing to Eq. fl25|), and recalling 
the expression for the Kraus operators, Eq. (||), it follows 
that this term is directly related to bath matrix elements 
of Hsb, which give rise to a Lamb shift 58 . When 
the system Hamiltonian is included, it appears in the 
term as well. However, recall that we are developing 
an approach that is explicitly model- independent. Hence 
the only quantities we will use are the QPT-measurablc 

Xai.O- 

We now wish to find an appropriate set of BB controls 
in order to eliminate the noise on our qubits. It should 
be clear from the discussion we just presented that this 
noise is unitary errors (and not decoherence) , since in the 
short-time limit relevant for BB we only deal with the 
bath-induced Lamb shift [decoherence arises from terms 
that are 0(t 2 )]. As noted above, from sections IIA and 
[IB we find that, under the action of BB controls, the K 



transform as 



Ka 



BB 



(30) 



(Here a denotes a vector of arbitrary index.) This trans- 
formation is the basis for much of what follows. 



C. Qubit Storage 

For the storage of information (without computation) 
in qubits, we need to preserve the density matrix under 
time evolution, so that p(t) — p(0). Let us denote BB- 
modified quantities by a tilde. In this case we should 
have, using Eq. (p9|), 



[5(t),p(0)]=0 



(31) 



as the BB control objective. Since S does not contain an 
identity component 1 we require that 



S(t) = 0. 



(32) 



We proceed to turn this into a condition on BB pulses. 

Recall that K$ 1 denotes an operator with exactly one 
non-identity term (one of the three Pauli matrices acting 
on an unspecified qubit). There are therefore 3N such 
operators, which we now denote explicitly as erf, where 
i = 1, N , and a — 1, 2, 3. Under the assumption of a 
linear system- bath coupling, Eq. (|24|), it is clear that the 
BB-pulses need only involve tensor products of single- 
qubit unitaries, i.e., 



N 



u k = < 



(k) 



Then Eq. @ becomes 

= pE^ fc)t ^ fc) . 



(33) 



At this point it is useful to again introduce real rotation 
matrices R to represent the BB-group: 



(34) 



Here i runs over qubit indices; k € {0, 1, \G\— 1}; R 1 '^ 
is in the adjoint representation of the group SU(2) [i.e., 
jj>i ; (fc) g 50(3)] acting on the i th qubit and has matrix 

elements R 1 '^ . Now let us consider the transformation 
of S(t) under the BB controls. To simplify notation let 
us denote 

C a = Im(xi, )- (35) 

Then from Eq. @: 

7S(*)«E£&tf = E?-^ ( 36 ) 
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where 



(fc) 



(37) 



(38) 



For storage we require £ = 0, i.e., 

$=Mj&o) = VAi- (39) 
Thus, solving for each i the set of linear equations 

EE Im (xL,o)< fe) =0 (40) 



/or t/ie rotation matrix elements R i '^ k \ in terms of mea- 
surable parameters Xaofthe output of a QPT experi- 
ment), determines the BB pulses empirically. The pulse 
form of the BB controls is determined through Eq. (34). 

Note that if lm(xa,o) = ^ m (Xa,o), i.e., there is no de- 
pendence on qubit index (collective decoherence f| B[), 



general conditions from the empirical BB perspective, 
which have the advantage that they can used to deter- 
mine the required set of BB pulses directly from a QPT 
measurement and a stipulated, wanted system Hamilto- 
nian. 

Let us consider the case of single-qubit operations first. 
In this case the system Hamiltonian need only contain 
a single non-identity operator (Pauli matrix) per qubit, 
as in Eq. (p4|). Therefore the development of the previ- 
ous subsection applies. The difference, however, is that 
now instead of the storage condition of Eq. (|3^) we re- 
quire the BB-modified x-matrix elements x% o t° assume 
values that correspond to a wanted evolution (or sys- 
tem Hamiltonian S w ). Let us denote the corresponding 
wanted (real)x-matrix elements by ufo (they can easily 
be calculated from a Hamiltonian - see below); then the 
empirical BB condition replacing Eq. (J40|) becomes: 



1 

W\ 



EE Im (xL 



,0) H a/3 



(42) 



This once again has to be solved for the rotation matrices 
Jl^(. k ) ) w ith elements R 1 ^ , given the empirical data x = 
Im(Xa o)- This too, can be written in a form amenable 
to a geometric interpretation: 



(fc) 



Im(x 



(43) 



then the same set of rotation matrices {R^}j?l (with 



matrix elements R^l) can be used for all qubits, as al- Now the average over the rotation matrices acts to rotate 



L ctf3 

ready pointed out in |]16| in terms of unitary BB controls. 
It also shows that, for complete symmetrization, one need 

only ensure that Ylk^aB ~ f° r au Q iA independent 
of the decoherence mechanism. 

Finally, note that we can rewrite Eq. (EH) as: 



the QPT output vector to a desired vector for the 
qubit, w l . Eq. ( |43| ) is our third main result. 



(41) 



In this manner it is clear that what we are looking for is 
a group of rotation matrices {R 1 '^}, acting on qubit i, 
whose average 

W\ Sfe Rl ' [k) acts t0 annihilate the QPT 

measurement output vector Im(x ) ■ This is the geometri- 
cal interpretation of the empirical BB condition. Eq. ( pd| ) 
is our second main result. 



D. Single-Qubit Operations 

Now suppose that we are interested in quantum com- 
putation. In this case we must allow for single- and 
two-qubit operations, such that these are not eliminated 
by the BB controls. In the model-dependent approach 
this translates into the (sufficient) condition that the BB 
generators commute with the Hamiltonian that is imple- 
menting the computation (l8| E9[. Here we derive more 



E. Two-Qubit Operations 

In order to implement two-qubit operations we must 
allow for a system Hamiltonian that contains two-body 
interactions. Therefore it is useful to comment on what 
happens when also the system-bath Hamiltonian contains 
higher order coupling, e.g., second order: 



N 



H 



(2) 
SB 



E ■ Gi - 



B, 



(44) 



j>i=l 



where Gy is a second-rank tensor. In this case both the 
third and fourth line of Eq. (|2^ ) contribute terms that are 
bilinear in the Pauli matrices, i.e., they contribute X)a 2 
and X^Sj 1 respectively. It is important to distinguish 
these bilinear terms from additional bilinear terms that 
arise when the expansion is taken to 0{t 2 ). The latter 

may arise from Hgg [Eq. (p4])j and will contribute to the 
non-unitary, decohering, part of the evolution. However, 
in the context of BB controls we are only interested in 
the ultra-short time limit 0{t), where the bilinear terms 
arising from H^only contribute to the Lamb shift. Thus, 
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the BB pulses that are appropriate for both two-qubit 
operations and a second-order system-bath Hamiltonian 
will be elements of SU(4). 

As before, the quantities extracted from the QPT mea- 
surements will be the imaginary part of the %-matrix 
which we abbreviate using a matrix £, as in Eq. (35). 
In this case the modified evolution will provide for the 
possibility of two-qubit interactions. Thus, generalizing 
from Eqs. (|9|),@: 



= -*i> ( 45 ) 



where 



(46) 



is a 4 x 4 matrix of coefficients. Under the action of the 
set of BB controls, 



\q\ 
E 

ij 



where 



k ij 



1 



\ "* r>ij;(.k) , pij 



(47) 



(48) 



and the rotation matrices R E 50(15) are defined 
through: 



U. 



(fc)t 



(a? ® a 



0\ Tj(k) = 
3 ) ij 



7(5 



Again, let us describe the target, or wanted, evolution 
by the ^-matrix w. In analogy with Eq. ( |34"| ) the matrix R 
is in the adjoint representation of the group and thus can 
be viewed as a rotation in the vector space of Hermitian 
matrices. In this case, only a subgroup of the rotation 
group 50(15) is represented by the adjoint action. (This 
is true for all SU(n), n > 3). The expression analogous 
to Eq. (Ea) becomes 



(k) Hj 



Wi 



or using explicit index notation, 



1 



— p % 3 T?3\V> 

\Q\ 

11 k -)8 



ij p*j;(fc) _ ...ij 



~ W af3- 



(50) 



(51) 



Thus the two qubit case involves solving for the 225 ele- 
ments of each of the rotation matrices R l ^''^ k \ given the 



QPT data £ y and the desired Hamiltonian Wij. After the 
rotation matrices are found, one obtains the BB pulses 
by inverting Eq. (pil) for the U 



(k) 



While this seems 



like a daunting task in general, it should be numerically 
trac table, and is illustrated for a simple example in sec- 
tion IV E| below. Eq. ( [sb"| ) is our fourth main result. 



F. Generalization to Encoded Qubits 

Before moving on to examples, we generalize the em- 
pirical BB condition to encoded qubits, such as arise in 
the theory of quantum error correcting codes (QECC) 
iPl 1^1 3; PI an d decoherence-free subspaces (DFS) || 0, g, 
PC flO[ |11|. In both cases it is highly desirable to let the 
experiment determine a "tailored encoding" , since the 
experiment knows the decoherence processes that govern 
the system we wish to protect better than any model one 
can design! Furthermore, combining QECC and DFS 
with the BB method has proven to be a powerful tool 
||, H H ||, §§. Now, both QECC and DFS can be de- 
scribed in terms of a stabilizer ^ [l0| . A stabilizer group 
for a set of codewords, i.e., a code space, is a subgroup 
(of the Pauli group for QECC, and of the of the group of 
all unitary transformations for DFS) that leaves the code 
space invariant. A code (whether QECC or DFS) can be 
completely specified in terms of its stabilizer |6Cj] . 

Let 5 be the vector space (group algebra) generated 
by real linear combinations of the set of generators of the 
stabilizer group. Any member of S will leave the code 
space invariant. Thus the most general "error" we can 
allow in the outcome of a BB procedure, when compared 
to a given wanted Hamiltonian, is so that the error is in S. 
Then the outcome of this "erroneous" evolution will be 
correct up to an overall phase. As before, let Im(xw) be 
the coordinates of the vector corresponding to the desired 
Hamiltonian evolution and Im(%) the actual vector after 
BB operations. Formally, the condition is: 



5 - S w — ih Im(x) - Im(jttii) ■ K e S 



(52) 



which should be compared to Eq. (|23|). This equation 
may be interpreted in one of two ways. First, given an 
encoding and a wanted Hamiltonian S w , it can be solved 
for the BB operations that are needed for the suppression 
of errors on the code subspace. Second, given a physically 
implementable set of BB operations, it can be solved for 
a compatible code (by finding the stabilizer) . Abstractly, 
this procedure may be seen as a projection of the open 
system evolution onto an evolution which is in the stabi- 
lizer of the code space. The geometric projection opera- 
tion completely reduces to the group-theoretical projec- 
tion onto the commutant given in [ 15[ i n the case that 
the set of RW , form a discrete group |p4| . The emphasis 
here is a geometric picture of the empirical operations 
projecting onto the stabilizer group of the code. Note 
that quite generally, Eq. (p2|) gives an empirical means 
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of identifying a subspace encoding such that the BB oper- 
ations drive the evolution into a subspace which does not 
affect the encoded states. This implies a general, empiri- 
cal means for the creation of a DFS |^7f . 

If the BB procedure is imperfect there will be an error 
component remaining. The error vector E is given by the 
difference between the BB-modified and wanted Hamil- 
tonians in the (n 2 — l)-dimensional vector space where 
our geometric picture holds: 



E = Im(x) - Im(xt, 



(53) 



The vector E gives the magnitude and direction of the 
error (i.e., the basis elements A, give the type of error, 
e.g., bit-flip and/or phase-flip, etc.). The corresponding 
scalar quantity is 

d(Imf,Im(x„) = [Tr(Im(f) - lm(x w )) 2 } 1/2 . 
The error (^3|) can be generalized to 

d(S, Im(x) ) = mm(B - Im(x) ), 

and likewise, the corresponding scalar quantity is 
d(S,lm(x) ) = mm[Tr(B - Im(|) ) 2 ] 1/2 , 

which can be visualized as in Figure |l|. 



(54) 
(55) 

(56) 




FIG. 1: Visualization of the error d that remains after the ap- 
plication of BB pulses. Im(x) is the BB-modified coordinate 
vector of the evolution, and d measures the distance to the 
closest element of the algebra of the stabilizer group of the 
code. 



Similar conclusions were presented for the unencoded 
case in [pi|. 



IV. EXAMPLES 

In this section we study a couple of simple examples 
that illustrate the formalism developed above. 



A. One-Qubit Example: Storing a Qubit in the 
Presence of Pure Dephasing 

Let us consider a simple model: a phase- flip error (pure 
dephasing) on a single qubit. To first order, this gives a 
density matrix of the form 



Ps 



igt 



(57) 



where the prime indicates the density matrix for the 
qubit after the interaction with a bath. The coupling 
constant g is a measure of the strength of the interac- 
tion. The bath time-scale is the inverse of the bath high- 
frequency cutoff, which is a separate parameter. Suppose 
that we wish to find a set of BB pulses that store this 
qubit. 

The first step in the empirical BB procedure is to 
measure the superoperator using QPT. Here we would 
discover that the interaction causes a phase-flip error 
which corresponds to K oc a z . I.e., a measurement of 
the Y-matrix would yield [by comparison of Eq. ( |57|) to 

Eq. ©]: {Im(xi,o)}<*=* )3 /,* = {°> °> -9/ 2 }- 

The next step is to find the optimal set of BB opera- 
tions. This we can do by solving Eq. ([II]) for the rotation 
matrices with the measured %-matrix. This yields 




IGl-l 

E 

k=0 



R 



(k) 

3/3 



0. 



(58) 



In accordance with the BB operations forming a discrete 
subgroup, k = corresponds to the identity. For n = I, 
since -R31 = for the identity rotation, 



+ 



I6I-1 

E 

k=l 



R 



(k) 

3/3 



|6| 



and I 



R 



33 



(59) 



fc=i 



where f3 = 1, 2. 

The best set of BB operations is the set that accom- 
plishes the task at hand and has the fewest elements \Q\. 
We now find a set with \Q\ = 2 (corresponding to a parity 
kick solution ||, H|). I.e. we seek a rotation matrix 



whose unspecified elements are arbitrary in as far as that 
they are not determined by the QPT data. To transform 
from the rotation matrices back to the BB pulses [i.e., 
from SO(3) back to SU(2)], we use the general result: 




a a cos(2#) + 2n a (n ■ (?) sin 2 (#) 

-{n x o%sin(20) (60) 



Here R s 50(3), n is a unit vector along the axis in K 3 
about which a rotation through an angle 9 is performed 
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(these 4 parameters parametrize the 50(3) rotation ma- 
trices). Let a = 3, then we know from the form of 
that: 

— a z — a z cos(2#) + 2n 3 (n ■ a) sin 2 (0) — (h x a) z sin(2#) 

It is simple to check that (mod27r) the unique solution to 
this equation is: 9 — ±7r/2, ti 3 = 0. This implies U — 
e ±m-crir/2^ w j^jj _ (ni,7i2,0), but otherwise arbitrary. 
The BB pulse thus must correspond to a rotation in the 
x — y plane on the Bloch sphere, which is the expected 
result as the error was along the z axis. 

It is likely that in a real experiment pure dephasing will 
not be the only source of decoherence. Let us consider 
a situation where this was the dominant source, so that 
our QPT measurement that yielded {ho-(Xa-o)}a=x,y,z = 
{0, 0, — g/2} actually contained an ^-component as well, 
which was too small to be noticed while the dephasing 
process was present, e.g., because the two errors may well 
have different characteristic time scales. Suppose that we 
perform another QPT measurement while applying the 
BB pulses found above (that eliminated dephasing) and 
find a residual error of the a x (bit-flip) type. This is an 
instance of a learning loop, which we discuss in Section 
below. 

In this case, consider the total Hamiltonian 

H = |-a x ® (I + a x ). (61) 

Proceeding in exactly the same manner as before we 
determine the required BB operations. We find that 
we need to implement U = e ±n <T7r / 2 j where now h = 
(0,712,713). Combining this and the condition n — 
(ni, tt.2,0), we find that we need to use h — (0,7i2,0). 
Thus bit and phase flips can be corrected using the cor- 
responding single BB operation, which is determined em- 
pirically from an experiment with a learning loop process. 
This is an optimal set since it will eliminate both errors 
with only one (non-identity) BB pulse per cycle. 

B. Two-Qubit Example: Computation Using the 
Heisenberg Interaction in the Presence of 
Independent Dephasing 

As indicated above, the problem in the two-qubit case 
can be quite involved since in general it requires finding 
the elements of rotation matrices in 50(15). To illustrate 
the formalism we consider a simple example. Suppose 
we wish to implement a Heisenberg exchange interaction 
J<i\ ■ 02 between the two qubits (Heisenberg exchange is 
important in a number of promising solid state proposals, 
and is an interaction that is all by itself universal for QC; 
see, e.g., |[l], and references therein). Then the wanted, 
Heisenberg interaction is determined from 

H Heis = Jai ■ a 2 = ^2 ^2 a ? ' w afi ' °f ' 

ij a (3 



so that it is described by the matrix 

/ 1 0\ 

W12 = J 1 . (62) 

V° 1/ 

Further, suppose that our QPT measurements suggest 
that the source of decoherence in the experiment is in- 
dependent dephasing on the two qubits. This will be 
detected through QPT by producing the following: 

*-S(t) « g x at + g 2 a* 2 =EE<" £p " 4 ( 63 ) 

ij a/3 

Independent dephasing will thus be described by the ma- 
trix elements 

£ 3 12 o=2i, Co 2 3=52- (64) 

To find the set of BB pulses we would now need to solve 
Eq. (|5|) for the rotation matrix elements, and then de- 
termine the corresponding SU (4) transformations, in a 
manner analogous to what we did above in the single- 
qubit example. As noted above (see also ^4|) solving 
these equations is not, in general, trivial. In this sim- 
ple example, however, a set of BB controls can be found 
noting the trace orthogonality of the two algebraic basis 
elements p4| , those corresponding to the exchange and 
those corresponding to the errors. Rather than going 
through a full derivation, we present the solution. To 
remove the independent dephasing through a parity-kick 
procedure, without affecting the Heisenberg exchange in- 
teraction, it is possible to use independent qubit interac- 
tions which form the following pulse 

U = UxU 2 = cxpHK + of )?r/2) = -o-fcrf. (65) 

By direct calculation one can show that 

[C/,H H cis] = and {C/,5} = 0. (66) 

The first commutation relation ensures that the parity- 
kick pulse can be applied during computation with 
-^Heis, while the second (anti-)commutation relation is 
the parity-kick condition [p[ [lj, |||. Thus the de- 
sired evolution is achieved. The pulse U is certainly not 
unique, and a general solution of Eq. (|m]) would yield a 
variety of other possible pulses. 

V. OPTIMIZATION ALGORITHMS 

As indicated in the single-qubit example discussed in 
the previous section, the empirical BB procedure can 
benefit from the incorporation of an off-line learning loop, 
that acts as an optimization algorithm for the BB pulses. 
Such learning loops have proven very successful, e.g., in 
quantum chemical applications, where they are typically 
used to optimize the yield of a chemical reaction, steer 
a system towards a desired state, or perform a cooling 
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task mm, n 00 ■ R ° u g hi y ; the idea is t° e uide 

a quantum system toward a desired goal by letting a 
learning algorithm optimize a classical control field (e.g., 
a laser pulse). An initial field is guessed and applied 
to the quantum system. The output is measured and 
input into a search algorithm (e.g., a genetic algorithm 
p7[ |62[ ) , which tries to optimize the field in order to get 
closer to the desired goal. The experiment is then re- 
peated with the new field, and the process is repeated 
until it converges to the desired goal to within a pre- 
scribed tolerance. 



A. Variational Optimization 

We first present an outline of a variational optimization 
procedure, which can in principle be used to tailor our BB 
pulses. Our presentation follows the standard approach 
in the quantum control literature, e.g., p5| , |36| . The 
general control problem can be stated as follows. We 
seek a system Hamiltonian, H Cl which modifies a given 
(total, system-bath) Hamiltonian H, so as to produce the 
desired effective Hamiltonian 



H(ui(t),t) and we may formally write the solution as 



H = H + H r 



(67) 



The control Hamiltonian H c may be composed of several 
possible terms, 



(68) 



where the ti»(t) are usually pulses in QC, analogous to the 
control fields in NMR and quantum optical systems. I.e., 
the Ui{t) are control fields that may be turned on and off 
as desired. The unitary evolution will proceed as usual 
according to U(t) = T exp(-i /* H(t')dt') with H, H c 
acting simultaneously. Thus controllability is determined 
by the group space that one is able to generate by the ex- 
ponentiated vector fields H 63 . For the robust storage 



of a qubit using BB controls we require the elimination of 
the interaction Hamiltonian H. This would correspond 
to having U « I . We also wish to use as few BB oper- 
ations as possible due to the time constraints. Thus we 
seek to minimize the difference between the BB-modified 
Hamiltonian S, and the desired Hamiltonian S w : 



(69) 



where S and S w are the appropriate modifications to Eq. 
(|l8|). E.g., for storage we would want S w = 0. One 
may now consider the standard controllability problem 
in terms of a desired state of the system, to be reached 
from some initial state \a Q ) = \a(t = 0)): 



i\a(t)) = H\a{t)) 



(70) 



(formally we should have included bath states as well, 
but we omit these for notational convenience). H = 



\a(t)) = U(t,to)\ao) 



Texp<j / H{ Ul {T),T)dT 

to 



M, (71) 



where T is the time-ordering operator (the H(ui(r), r) do 
not necessarily commute). This can be seen as essentially 
a Heisenberg picture control problem and one can 
thus eliminate the direct inclusion of the state itself. The 
short-time approximation enables us to remove the time- 
ordering and write the expanded form 

U(t,t ) = lim [exp{-iH N -it N - 1 }...exp{-iH Q t }] 

(72) 

At this point we may invoke the assumptions of the BB 
operations that they be short, strong pulses and the evo- 
lution in between them be that of the free system-bath. 

From this one may also see the connection with the 
standard control theory that often uses the final state as 
the "output" of the control. This may be used for nu- 
merical algorithms which are associated with a learning 
or realtime feedback loop. However, it is clearly desir- 
able to have both the operator and state pictures (i.e., 
with and without the explicit state dependence), since 
one may often wish to consider the control of the evolu- 
tion rather than that of the state. In fact, in quantum 
computing, the control objective is noiseless evolution 
rather than simply obtaining a target state. 

To optimize the BB procedure, the difference between 
the BB modified controls and the ideal evolution should 
be minimized. This may be achieved in the continuum by 
solving the variational problem with a variable end point. 
The appropriate variational problem can be formulated 
as the minimization of a cost function J (3^, [36l |3^], |3^. 
fl(i| , expressed in terms of the control fields{tti} and 
cycle time T c as: 



J = 



MT C 



(TT{[S(u(r),T)-S w (r)} 2 }) 1/2 dT, (73) 



where a 2 = a^a and we have used M cycle times T c for 
the end point (which is not fixed). One may add exper- 
imental constraints, such as finite pulse energy, smooth- 
ness of the pulse shapes, etc. |35|, [56|. This is a standard 
variational problem for which we would seek S J = and 
5 2 J < 0. The outcome, i.e., the solution to the varia- 
tional problem, will be the optimal control fields 
Note that these fields will be approximately continuous 
for large M and small At. Then we may approximate 
them by a discrete set of BB operations (traditionally dc- 
fined as piecewise continuous controls; see, e.g., HJQ). 
However, it is to be expected that one of the advantages 
of the optimization procedure is that it will yield pulses 
that are easier to implement physically than the pulses 
coming out of a standard BB analysis, since the optimiza- 
tion procedure can be formulated to explicitly take into 
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FIG. 2: The learning loop control diagram for TQEC. 



account experimental constraints. In fact, experience 
in quantum chemistry shows that the pulses found by 
an optimization procedure are often highly non-intuitive 

Questions of convergence, etc., can be avoided by the 
use of a small number of cycle times. This will reduce the 
problem, under the BB assumptions, to a search on a dis- 
crete space. This is the space of discrete, or finite order, 
subgroups of unitary groups. Fortunately, for quantum 
computation, we require only one- and two-qubit oper- 
ations which reduces our search spaces to those of the 
discrete subgroups of SU{2) and/or SU(4). These have 
recently been classified (see |66[ and references 
therein). We will not pursue the variational formulation 
further here. An actual variational optimization calcula- 
tion will be presented in a future publication. 



B. Learning Algorithm 

In certain cases it may be possible to perform a large 
number of experiments on identically prepared samples, 
differing in the applied control fields. In this case, in- 
stead of solving a variational problem to find optimal BB 
pulses, one can try to let the experiment guide an off-line 
learning algorithm (typically a genetic algorithm) to an 
optimal solution |37| . This algorithm is part of a learning 
loop, described in Fig. |[ 



The learning loop consists of the following steps, which 
are repeated iteratively in the learning process. 

1. A quantum state is input for a particular informa- 
tion processing task. 

2. The state is allowed to interact with a bath and 
undergo noisy evolution in the experiment Here 
we may choose to apply BB pulses to modify the 
evolution. 

3. The resulting evolution is obtained through quan- 
tum process tomography, QPT. 



4. The QPT data is analyzed by the learning algo- 
rithm to find an improved BB strategy. This in- 
volves solving the key Eqs. (p3|) , (|50|) . 

5. The previous steps are repeated until convergence 
to within a prescribed tolerance. 

The result of the procedure is an optimized set of BB 
pulses. This set includes (i) the least number of BB op- 
erations that will reduce or eliminate the noise in the 
system, and (ii) the optimal ordering for this minimal 
set. 

Let us emphasize that: (1) no knowledge of the total 
Hamiltonian or noise process is assumed (i.e., the opti- 
mal implementation is determined empirically), and (2) 
no assumption is made about the quality of the BB op- 
erations, only that they should improve the fidelity of 
the desired operations. (Of course we know from earlier 
work H H 0, P, [if 0, H, H that the BB oper- 
ations should be implemented as strong fast pulses, but 
imperfect implementation will still reduce noise.) 

Finally, let us note that the learning process could in 
princi ple be incorporated in a real-time feedback loop 
(e.g., Q and references therein), but this would require 
a very fast numerical algorithm to solve Eqs. (f43|),(50). 



VI. CONCLUDING REMARKS 

In order for methods that reduce decoherence and noise 
in quantum information processing tasks to succeed in 
the real world, they must be confronted with experi- 
mental data, and allowed to be optimized in response 
to this data. This is the approach we have taken here, in 
the context of the dynamical symmetrization, or "bang- 
bang" (BB) method. We have developed a formulation of 
the BB method that allows one to tailor the BB control 
pulses in response to data acquired by a quantum process 
tomography experiment. The experiment supplies a set 
of numbers that characterize the noise processes occur- 
ring on a short time-scale. From these numbers one can 
determine an optimal set of BB pulses, by solving a set 
of linear equations, in particular Eqs. (p3|) ,(41), (p3|) , (^Oj) - 
These equations correspond to different tasks one may 
wish to implement with the help of the BB pulses (re- 
spectively, general storage, single-qubit storage, single- 
computation, two-qubit computation), and yield a set of 
rotation matrices (R) that correspond to BB pulses that 
perform the desired tasks. 

A promising generalization of a single-shot 
tomography-BB experiment is to introduce an off- 
line learning loop, that uses the above equations in 
order to determine an optimized set of BB pulses. The 
learning process incorporates tomography measurements 
from a previous round in order to find improved BB 
pulses for the next round. We have briefly discussed how 
such a loop, and a concomitant variational optimization 
procedure, can be designed. 
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Throughout this work we have emphasized that our 
results have an intuitive interpretation in terms of a ge- 
ometric picture, wherein the effect of BB pulses is to 
rotate a coordinate vector representing a noisy Hamil- 
tonian to a desired Hamiltonian. The geometric picture, 
via Eq. (|52|) , also enables the determination of the ability 
to create an encoding (such as a decoherence-free sub- 
space) using empirical data and the available set of BB 
pulses. Alternatively it can be used to describe the ap- 
propriate set of BB operations required to eliminate noise 
from an encoded set of qubits. 

We hope that the results presented here will stimulate 
experiments in which real data will drive the determina- 
tion and application of appropriately tailored BB pulses. 
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